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We study the properties of strange quark matter in equilibrium with normal nuclear matter. 
Instead of using the conventional bag model in quark sector, we achieve the confinement by a 
density-dependent quark mass derived from in-medium chiral condensates. In nuclear matter, we 
adopt the equation of state from the Brueckner-Bethe-Goldstone approach with three-body forces. 
It is found that the mixed phase can occur, for a reasonable confinement parameter, near the normal 
nuclear saturation density, and goes over into pure quark matter at about 5 times the saturation. 
The onset of mixed and quark phases is compatible with the observed class of low-mass neutron 
stars, but it hinders the occurrence of kaon condensation. 

PACS numbers: 21.65.Qr, 26.60.-c, 25.75.Nq, 26.60.Kp 
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I. INTRODUCTION 

By far the study of neutron stars has been mainly fo- 
cused on the relationship between the equation of state 
(EoS) of nuclear matter and the observed maximum 
mass. The connection has been achieved by solving the 
hydrostatic equilibrium equations based on general rel- 
ativity. The first generation observed masses exhibited 
an average value around 1.5 solar masses. This value re- 
quires a soft EoS that can be easily obtained by introduc- 
ing new degrees of freedom like hyperons, kaons or quarks 
accompanied or not by a phase transition. Sometimes 
the softening was so large that the neutron star is pre- 
dicted to collapse into a black hole, as for the SN1987A 
[l| . In the new generation of observations the masses are 
distributed within a large range, up to 2 solar masses, 
that requires a stiff EoS, i.e., hadronic matter without 
new degrees of freedom. Since it is hard to imagine pure 
hadronic matter to sustain the high pressure predicted in 
the inner core, new scenarios have to be advanced to ex- 
plain the coexistence, in the phenomenology of neutron 
stars, of low and high mass spectra. 

Recently [2| it has been argued that the two observed 
classes of neutron stars might correspond to two different 
evolutionary scenarios of neutron stars. In one case, the 
hot and dense remnant of the supernova explosion rapidly 
evolves into a hybrid star, where the transition to a quark 
phase softens the nuclear matter so that M « 1.5M Q ; in 
the other case a slow evolution could lead the neutron 
star to a large mass via a mass accreting from the cou- 
pling with a white dwarf. From this point of view the 



destiny of the remnant is strongly affected by the ini- 
tial conditions, i.e. density, temperature, leptonization 
degree etc. For instance, if the mass of the remnant is 
below the mass threshold for quark nucleation the tran- 
sition to the quark phase is forbidden Q. If the mass 
is slowly accreting the transition is allowed. The role of 
temperature or other parameters defining the initial state 
of a new born neutron star has not yet been studied. 

To investigate the possible phase transition to quark 
matter in neutron stars, we need also to know the EoS of 
quark matter. Although we have in hand the fundamen- 
tal theory of strong interactions, i.e., Quantum Chromo- 
dynamics (QCD), we still do not know the true ground 
state. It is now generally expected that quark matter is 
in the color-flavor locked phase (CFL) [J] at extremely 
high densities when the finite current mass of strange 
quarks becomes unimportant. In the density range from 
nuclear saturation to CFL, there may exist a rich and var- 
ied landscape of phases, e.g., the 2SC, g2SC, gCFL etc. 
Presently, however, these phases suffer from the so-called 
chromomagnetic instability problem for both the two- [5j 
and three-flavor [a, S @| cases. On the other hand, ex- 
periments show that quarks become asymptotically free 
rather slowly [9(. Therefore, in the present study we are 
deal ing with the ordinary strange quark matter (SQM) 

EEL 

The special problem in studying the EoS of ordinary 
quark matter is to treat quark confinement in a proper 
way. In the conventional standard approach, an extra 
constant term, the famous bag constant B, is added to 
the energy density of the system, which provides a neg- 



ative pressure to confine quarks within a finite volume, 
usually called a 'bag'. The quark mass is infinitely large 
outside the bag, and a finite constant within the bag. A 
vast quantity of investigation have been performed within 
the framework of the bag model [12| • 

As is well known, however, particle masses vary with 
environment. Such masses are usually called effective 
masses. Effective masses of hadrons and quarks have 
been extensively discussed, e.g., within the Nambu-Jona- 
Lasinio model [f3[ and within a quasi-particle model [lj] . 
In principle, the density dependence of quark masses 
should be connected to the in-medium chiral condensates 

mill- 

Taking advantage of the density dependence, one can 
describe quark confinement without using the bag con- 
stant. Instead, the quark confinement is achieved by the 
density dependence of the quark masses derived from in- 
medium chiral condensates [lj, [l8|, [lj|. The two most 
important aspects in this model are the quark mass scal- 
ing [17], Qj| an d the thermodynamic treatment [lg, [lj| . 
Both aspects will be reviewed in this paper. 

In the present contribution, the transition from hadron 
phase (HP) to strange quark phase (SQP) in the inner 
core of a neutron star is investigated within the fully 
consistent nuclear and quark models. In the hadron 
sector we adopt the equation of state from Brueckner- 
Bethe-Goldstone (BBG) approach with three-body forces 
(TBF) (HIHIll. This theory, being a completely mi- 
croscopic approach, can easily incorporate degrees of free- 
dom such as nucleon resonances [A(1232) or AT* (1440)], 
which are expected to appear at higher hadron densities. 
It is found that the mixed hadron-quark phase can occur, 
for reasonable values of the confinement parameter, a lit- 
tle above the normal saturation density, and can undergo 
the transition to pure quark matter at about 5-6 times 
the saturation. This result is quite different from the pre- 
vious results from Nambu-Jona-Lasinio (NJL) model in 
which the mixed quark phase can not appear at neutron- 
star densities [23 [2J]. Afterwards, the influence of the 
mixed and quark phases on the structure of compact stars 
is discussed by solving the Tolman-Oppenheimer-Volkov 
(TOV) equation and extracting the mass-radius plots for 
neutron stars. Finally, it is shown that the transition to 
the deconfined phase turns out to be incompatible with 
the onset of kaon condensation. 



II. EOS OF QUARK MATTER 

SQM has been one of the hot topics in nuclear physics 
since Witten's famous stability conjecture [ljj]. In many 
studies, the quark confinement was treated adopting the 
bag mechanism [ll|, [25J . An alternative approach to ob- 
tain confinement is based on the density dependence of 
quark masses 26]. This mechanism has been exten- 
sively applied to investigating the properties of SQM 
[27ll2a .l29l.l30j]. In this section, we first give a short review 
of the two most important aspects, with the main incon- 



sistency of the original model. Then we present a fully 
self-consistent thermodynamic treatment. The proper- 
ties of SQM will be given in the new treatment. In the 
present paper, however, the main application of the new 
approach is the study of the phase transition in compact 
stars. 



A. Confinement by density-dependent masses 

As mentioned above, the quark confinement in this 
model is achieved by the density dependence of quark 
masses. Therefore, the first important question is how 
to determine the quark mass scaling which can reason- 
ably produce confinement. Originally, the interaction 
part of the quark masses was assumed to be inversely 
proportional to the density [26|, [27j . This linear scal- 
ing has been extensively ap plie d to studying the prop- 
erties of SQMfplpl IM HJ- There are also other 
mass scalings [3l|, [32 1 . Their main drawback is that they 
are pure parametrizations without any convincing deriva- 
tion. Therefore, a cubic scaling was derived based on the 
in-medium chiral condensates and linear confinement at 
both zero [171 ] and finite temperature J19l] . This new scal- 
ing has been applied to investigate the viscosity of SQM 
and the damping time scale that is due to the coupling of 
the viscosity and r mode [331 ] , the quark-diquark equation 
of state and compact star structure [34], the properties 
of strangelets versus the electric charge and strangeness 
[33, and the new solutions for CFL slets yo|. In the 



present paper, we use the chirally determined quark mass 
scaling (l7l fL9l ] to study the phase transition in neutron 
stars. For this we need a completely self-consistent ther- 
modynamic treatment of the EOS of quark matter. 

The thermodynamic treatment of the system with con- 
finement via the density dependent quark masses has 
been a long story. Originally, the thermodynamic for- 
mulism was regarded the same as the constant-mass case 
[27J. In this first treatment, the internal pressure can 
not be zero, and the properties of SQM were rather dif- 
ferent from those in the bag model. But it was later 
pointed out that the difference was caused by the incor- 
rect thermodynamic treatment [281 ] . It was found that an 
additional term is to be added to the pressure and energy 
expressions [28J . This second treatment makes it possible 
SQM to be self-bound. However, two serious problems 
came out: one is the unreasonable vacuum limits, the 
other one is the discrepancy between the energy mini- 
mum and zero pressure. It was shown that the added 
term in the pressure, due to the density dependence of 
quark masses, should not be appended to the energy. Af- 
ter discarding this term in the energy, while keeping it 
in the pressure, the two inconsistencies mentioned above 
were immediately removed [181 ] . This third treatment has 
recently been extended to finite temperature (l9| . The 
thermodynamic formulism in Ref. [181 ] was also adopted 
in Ref. [30j, though a different quark mass scaling was 
used there. 



A common feature of the last two thermodynamic 
treatments [111 Il9l. 12811 , as well as other recent references 
using this model [30] , is that they all regard the thermo- 
dynamic potential the same as in the Fermi gas model. 
Because of the additional term, the pressure becomes ob- 
viously not equal to the minus thermodynamic poten- 
tial density, contradicting the thermodynamic equality 
P = — il for a homogeneous system. One can also eas- 
ily check that the fundamental differentiation equality 
(IE = J2i /^idrij for homogeneous systems at zero tem- 
perature was not fulfilled in the mentioned references. 

In the rest of this section, we will present a fully self- 
consistent thermodynamic treatment of the confinement 
by density-dependent mass model (CDDM). 



In the following, we show that the density-dependence 
of particle masses will modify the Fermi momentum, i.e. 
the relation in Eq. $5§ for free-particle systems should be 
modified to include interactions. In fact for the quark 
flavor i we have 
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Since the quark masses are density dependent, the deriva- 
tives generate an additional term with respect to the free 
Fermi gas model. We get 
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(6) 



B. Self-consistent thermodynamics in CDDM 

Let us consider a quark model with three flavors. De- 
noting the Fermi momentum in the phase space by r/j, 
the particle number densities can then be expressed as 



rii = gt 



(2-Khf 2tt 2 J 
and the energy density as 



d 3 P 9i I " 2 - 9i"i , , . 

= TT3 \ V dp= -^±, (1) 
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E = Y,§2J^'\/p 2 + m *P 2d P- ( 2 ) 



Equations fl) and (2]| are familiar expressions, where 
the summation index goes over all considered particle 
types. To let the model be valid for both particles and 
anti-particles, the particle number density, or accord- 
ingly, the Fermi momentum, are formally assumed to be 
negative for anti-particles. Therefore, in the upper limit 
of the integration, the absolute value has to be taken. 

If the particle masses rrii are constant, the relation be- 
tween the Fermi momenta Vi and the chemical potentials 
Hi is 



M? - wf or m = 



(3) 



As is well-known, however, the quark mass depends on 
density and temperature. In principle, the quark mass 
scaling should be determined from QCD, which is ob- 
viously impossible presently. Based on the in-medium 
chiral condensates, a cubic scaling law was derived at 
zero temperature [171 ] , and it has been recently extended 
to finite temperature [l9J . At zero temperature, we have 
the simple cubic scaling 
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where m q o is the quark current mass, n is the total baryon 
number density, the exponent of density is z — 1/3 [l7( |. 
and the constant D is to be discussed a bit later. 



where 
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with J7o being the free-particle contribution: 
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Due to the additional term in the chemical potential, 
the pressure also has an extra term. The inclusion of 
such a term guarantees that the Hughenoltz-Van Hove 
theorem is fulfilled in the calculations. 

In the quasi-particle model [371 ] . one also has an ex- 
tra term. Because the quark masses there depend on 
chemical potentials, and the extra term is not used in 
the relation between the Fermi momenta and chemical 
potentials, an effective bag constant has to be added to 
the energy expression to consider confinement [T3 |. 

In CDDM quark model, however, we do not need a bag 
constant any more. Quark confinement is achieved au- 
tomatically by the density dependence of quark masses, 
or by the strong interaction between quarks. In fact, the 
exponent z — 1/3 in E q. (f4l) is derived from the linear 
confinement interaction [171 ] . 

In the present model, the parameters are: the electron 
mass m e = 0.511 MeV, the quark current masses m U Q, 
m d0i m sCb an d the confinement parameter D. Although 
the light-quark masses are not without controversy and 
remain under active investigations, they are anyway very 
small, and so we simply take m U Q = rrido = 0. The current 
mass of strange quarks is 95 ± 25 MeV according to the 
latest version of the Particle Data Group (38| . 
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FIG. 1: Confinement constant range determined by stability 
arguments. Only in the region bounded by the full and dashed 
lines, SQM is absolutely stable. 



Conventionally, the stability of strange quark matter 
(SQM) is judged by the minimum energy per baryon [101 , 
IllllHQ- K it is less than 930 MeV (the mass of 56 Fe 
divided by 56), then SQM is absolutely stable. If it is 
bigger than 930 MeV but less than 939 MeV (the mass 
of nucleons), then it is metastable. Otherwise, if it is 
larger than 939 MeV, SQM is un-stable. In case of two- 
flavor quark matter, it should be no less than 930 MeV, 
in order not to contradict standard nuclear physics. This 
is the Witten-Bodmer hypothesis [l(J, [ijj • 

in Fig.[TJ we show the different regimes in the \^D-m s o 
plane. The area below the full line is forbidden where the 
energy per baryon of two-flavor quark matter is less than 
930 MeV. Above the dotted line, the energy per baryon 
of SQM is more than 939 MeV, and thus SQM is unsta- 
ble. The area bounded by the dotted and dashed lines is 
the metastable region where the energy per baryon is be- 
tween 930 MeV and 939 MeV. Only when the (£> 1/2 , m s0 ) 
pair is in the range between the full and dashed lines, 
SQM can be absolutely stable, i.e., its energy per baryon 
is less than 930 MeV. Therefore, the range of D values 
is very narrow for a chosen m s o value, if the Witten- 
Bodmer hypothesis is correct. If we take the modest 
value m s0 = 100 MeV, for example, then D 1 / 2 is in the 
range of 158-160 MeV. The lower bound 158 MeV is ob- 
tained by taking m U Q = m d Q = 0. If m u o and m d Q are 
given a small finite value, the lower bound can then be a 
little bit smaller, e.g. 156 MeV [HI]. 

Unfortunately we presently do not have a definite con- 
clusion on the stability of SQM, so we treat Dasa free 
parameter. However, the first condition, i.e., D greater 
than about (158 MeV) 2 , should always be satisfied. On 
the other hand, we can connect D to the pion mass m^, 
pion decay constant f^, pion-nucleon sigma term ctn, 
string tension ao, and the vacuum chiral condensate (qq)o 



D = 



3(2/ 7 r) 1 /V m 2 / 2 



(10) 



From the known range of the vacuum condensate, we 
can have an upper bound (270 MeV) 2 . Therefore, D 1 / 2 
should not be out of the range (156, 270) MeV. 



C. Properties of strange quark matter 

As usually done, we consider SQM as a mixture of u, d, 
s quarks, and electrons. The relevant chemical potentials 
Mu> Md, Ms; an d Me satisfy the weak-equilibrium condition 



Mu + Me = Md) 

Md = Ms- 



(11) 
(12) 



Because all particle masses do not depend on the den- 
sity of electrons, i.e., drrij/dn e = 0, Eq. ^ gives 



Hi = y(TT 2 ni) 2 / 3 + m 2 - m 



with 



j—u,d,s J 



for i = u, d, s quarks, and 

Me = V (37T 2 n e ) 2/3 + TO 2 



(13) 



(14) 



(15) 



for electrons. 

In Eq. (|14p . mi is the second term on the right hand 
side of Eq. ((4|, so we have dmi/dn = —zD/n z+1 = 
—zmi/n. The pressure is then obtained from Eq. (JSJ) as 



P = -Cl + n h - — S~] rij /( 

dn b . ^— ' V m , 

j—u,a.s 



(16) 



Substituting these expressions into Eqs. (fTTj) and (fT2"|) . 
we have 



^J(ir 2 n u f 3 + ml + ^J(3ir 2 n e ) 2/3 + m 2 

= vV^) 2/3 +m 2 . (17) 



and 
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We also have the baryon number density 



1 



n = o( n « + n d + n s ) 



and the charge density 



1 1 



! q = ^n u - -n d - -n s - n e 



(18) 



(19) 



(20) 



The charge-neutrality condition requires Q q = 0. 

For a given total baryon number density n, we can ob- 
tain the respective n u , rid, Ti s , and n e by solving the four 
equations (fTTj). ((I5]l. (fT9|) , pjf . The chemical potentials 
Mjo Md, Ms, and /i e can then be calculated by Eqs. (fT"3|) 
and (JTSJ). Therefore, the energy density of the quark mat- 
ter is a function of the baryon number density n and the 
charge density Q q , i.e., E q — E q (n, Q q ) or 



dE q 
dn 



dE* 
dQ c 



dQ 



q- 



(21) 



where the two partial derivatives, dE q /dn and dE q /dQ q 
are called the baryon chemical potential and charge 
chemical potential, respectively. It can be easily shown 
that they are connected to the quark chemical potentials 

by 



dE q 
dn 



fJ u 



2\JL a 



dE q 
dQ q 



= Hu- Vd- 



(22) 



In fact, according to the fundamental differentiation 
equality of thermodynamics, we have 



dE q = fi u dn u + fj, d dn d + n s dn s + fi e dn e 



(23) 



On the other hand, we have \x e = /J>d ~ M«; Ms = Mrf? 
n s = 3n — n u — nd, and n e — n u — n — Q q from Eqs. pip . 
([T2"]) . dn]), and ((201) ■ Substituting these four equalities 
into Eq. ([23)) leads to d-E = (/i u + 2^£;)dn +(//„ — Hd)dQ q . 
Comparison of this with Eq. (f2"T| immediately gives 
Eq. " 
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FIG. 2: Quark fraction vs baryon number density for D ' 2 = 
160 MeV and m s o = 80 MeV. 

In Fig. [51 the quark fractions, i.e., n„/(3n), nd/(Zn), 
ris/(3n), and the 10 times the electron number divided 
by the total quark number, 10000n e /(3n), have been 



shown versus the baryon number density for D 1 ' 2 = 160 
MeV and m S Q = 80 MeV. It is seen that the fraction of 
up quarks is nearly always one third. The fraction of 
down quarks increases rapidly with decreasing densities, 
while the fraction of strange quarks approaches to zero 
when the density decreases to a certain lower density. 
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FIG. 3: Comparison between the Fermi momenta and chem- 
ical potentials. Parameters are the same as for Fig. [2] 

In order to compare the relation between the Fermi 
momentum and chemical potential, we plot, in Fig. [31 
the Fermi momentum of up (solid line), down (dotted 
line), and strange (dashed line) quarks, respectively, as a 
function of the corresponding quark chemical potential, 
in both the present model (lines with a solid circle) and 
the previous model (lines with an open circle) . It is very 
obvious that the difference is very large, especially at 
comparatively lower densities. In both models, the Fermi 
momentum of up or down quarks is higher than that of 
strange quarks due to the fact that strange quarks are 
heavier than up or down quarks. For the same chemical 
potential, however, the Fermi momentum in the present 
model is generally bigger than that in the previous model, 
due to the quark mass-density-dependence which reflects 
the strong interaction between quarks. 

Figure [3] shows the energy per baryon of quark matter 
for different parameter sets in the present model. Each 
line has a minimum, corresponding to the lowest energy 
state (marked with a solid triangle). One can see that 
the pressure at this minimum is exactly zero. So this 
special point is marked with an open circle as well. At 
the same time, we also display the energy per baryon for 
the two-flavor quark matter by a dash-dot-dot line. We 
see that the two-flavor quark matter is less stable than 
SQM. Even for a smaller D value, e.g., (160 MeV) 2 , its 
energy will finally exceed that of SQM for a bigger D 
with increasing densities. 
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FIG. 4: Energy per baryon of quark matter in the present 
model. The parameter pair (D 1 ' 2 ,m s o) 



in MeV for the 
solid, dashed, dotted, and dash-dotted cures are (200,150), 
(180,120), (170,95), and (160,80), respectively. The dash-dot- 
dot line is for the two-flavor quark matter at D 1 ^ 2 = 160 MeV. 
It is very obviously shown that the energy minimum, marked 
with a full triangle on each line, is located exactly at the same 
point of the zero pressure indicated by an open circle. 



In principle, the CDDM quark model contains more 
physics than the simple bag model. To demonstrate 
this we plot, in Fig. [51 the baryon chemical potential 
in CDDM with D 1 ' 2 = 180 MeV and m s0 = 120 MeV, 
and in the bag model with B 1 / 4 = 180 MeV. In CDDM, 
the baryon chemical potential decreases with decreasing 
density to a certain value depending on D and m s o, then 
it increases very rapidly, i.e., it saturates at a definite 
density marked with a bullet. When the density is lower 
than the bullet, the derivative d 2 E/dn 2 becomes nega- 
tive, and so quark matter is unstable against phase sep- 
aration and falls apart at lower densities. In the bag 
model, however, the baryon chemical potential is always 
a monotonic function of density, which means that quark 
matter does not fall apart at any lower densities. The 
velocity of sound has also been plotted on the right axis. 
We observe that it is, in the bag model, nearly the same 
as that of a non-interacting Fermi gas. But in CDDM, 
it decreases to zero with decreasing densities. At high 
densities, it becomes asymptotically identical to the ul- 
trarelativistic limit, as expected. 



D. Quark matter at finite temperature 

Since the single particle energies depend on density and 
temperature via the quark masses, the thermal proper- 
ties should be founded on the canonical ensemble, but, as 
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FIG. 5: The baryon chemical potential in CDDM (VD = 
180 MeV, m a0 = 120 MeV) and in the bag model (B 1/4 = 
180 MeV). The former has a minimum at a lower density 
depending on the value of D and m s o, while the latter is 
always a monotonic function of density. The velocity of sound 
in both models has also been given on the right axis. 



is well known, the partition function is not easy to calcu- 
late. Therefore a different statistical procedure is usually 
adopted, which is based on the quasiparticle assumption. 
According to that the energy density is written as 



£ = E^Evp 2 + m ' f *feT), 



where the Fermi distribution function is 
fi(p,T)= ' 



l + e h(p,T)- M ,]/T- 



(24) 



(25) 



If antiparticles are included, the sum must be extended 
to antiparticles for which fii must be replaced by — \ii. 
From the Landau definition of the single-particle energy 
extended to finite temperature, we have 
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£i(p) - Mi, 



(26) 



where £i(p) = \Jp 2 + m 2 is the dispersion relation of free 
particles. The extra term fj>i can be added to the chemical 
potential, so defining 



m + Mi- 



(27) 



Accordingly, the net density of the particle type i is rii = 
9i E P IHp, t ) ~ k{P, T )} , or, explicitly, we have 



9i 



1 



1 



1 + e N(p)-M*]/r 



p 2 dp 
1 + e N(p)+/4]/T / 2tt 2 ' 
(28) 



Inverting this equation, one determines /z* as a function 
of n, so that the free energy density 



F = Y / F t (T, f ilm l 



E 



F ; H 



f: 



(29) 



with 



Ft = 9i 



Tin 



± 



e -(\/P 2 + m ?TM*)/T 

p 2 dp 



/; 



l + e (Vp 2+m ?TK)/ T 



2tt 2 



(30) 



will be a function of respective particle densities, instead 
of chemical potentials. One can then determine the real 
chemical potentials and pressure, according to the well- 
known relations 

dF 



Hi = 



dn^ 



P=~F 



E 



Him. 



(31) 



These quantities will completely describe the thermal 
equilibrium of pure quark phase and the transition to 
the quark-hadron mixed phase. A more detailed anal- 
ysis, where the thermodynamic formulism is developed 
to much more extent, is reported in Appendix [AJ In 
Ref. [19l ]. the j«j should be implicitly understood as fj,*. 




d ENsity 



20 



FIG. 6: Density and temperature dependence of the energy 
per baryon for the parameters \AD = 180 MeV and m s o — 120 
MeV. 

In Fig. [SI we give the energy per baryon of SQM as a 
function of both density and temperature. The parame- 
ters used for this three-dimensional plot are \f~D = 180 
MeV and m s0 = 120 MeV. 



III. NUCLEAR MATTER IN BRUECKNER 
THEORY WITH THREE-BODY FORCES 

The Brueckner-Bethe-Goldstone (BBG) theory is 
among the most advanced microscopic theories of nu- 
clear matter. In the recent years it was recognized that 



the three-body forces, which are expected to have a dom- 
inant role at high nuclear density, also affect the satura- 
tion point and, in fact, after including the three-body 
forces in the Brueckner theory, the em piri cal saturation 
properties are reproduced quite well [2l], (39|. In the 
neutron star interior, where high baryonic density val- 
ues are reached, processes like the excitation of nucleon- 
antinucleon pairs (Z-diagrams) and nucleonic resonances 
(together with the production of other hadrons) sizeably 
influence the two-body nuclear interaction. The former 
process involves the virtual excitation of negative energy 
states, which is absent from the standard Brueckner the- 
ory, and thus it represents a pure relativistic effect. From 
the comparison with the Dirac-Brueckner theory it turns 
out that it is by far the most important relativistic ef- 
fect [33. This process together with nucleon resonances 
can be incorporated in the interaction as medium virtual 
excitations in TBFs. One could guess that many body 
(more that three) forces are also important at high den- 
sity as large as p sw 1 fm _ , but it is hard to imagine that 
pure baryon matter can exist at so high density. At lower 
density two and three-body forces are dominant since the 
hole line expansion is, roughly speaking, an expansion in 
density powers. 

The global effect of TBFs at high density is strongly 
repulsive, leading to a remarkable increase of the max- 
imum mass in the study of the neutron star structure. 
But also a correct estimate of the saturation point is 
important, since, as we will see below, in strongly asym- 
metric nuclear matter the threshold for the transition to 
mixed nucleon-quark phase can appear close to the satu- 
ration density. Therefore the corresponding EoS could be 
used as input for transport-model simulations of heavy 
ion collisions, where strongly isospin-asymmetric systems 
are formed in central events. 



A. BBG equations 

The Brueckner theory extended to TBF is described 
elsewhere [2(J, |2l|. Here we simply give a brief review 
of the Brueckner-Hartree-Fock (BHF) approximation at 
finite temperature T [4CJ, [4j], [42[ . The starting point is 
the reaction G-matrix, which satisfies the BBG equation, 



n , ^ V" |fcifc 2 )Qfc 1 ,fc 2 (T)(fcifc 2 | ^ „s 

G{uj,T) =i; NN + u NN > — G{u,T), 

fcfa ^-e fcl (T)-e fc2 (T) 

(32) 
where fcj = (fci,0i,Tj), denotes the single particle mo- 
mentum, the ^-component of spin and isospin, respec- 
tively, and to is the starting energy. The G-matrix, 
the Pauli operator Q and the single particle energies 
efe(T) = k 2 /2m + Uk{T) depend on the neutron and 
proton densities and temperature. The interaction dnn 
given by 



wnn 



ybare 



VI 



c rr 



(33) 



where V 2 barc is the bare two body force (2BF) and V£ s is 
an effective 2BF derived by the average of the bare TBF 
on the third particle as follows 



llr^ /drVr 3 'C(r 3 ')[l - ^(r' 13 , T)][l - V (r' 23l T)] 

n 

xW^P^\Pir^)(j> n {n)[l - J7(ri3,r)][l 



v(r23,T)]. 

(34) 

Since the defect function r)(r, T) is directly determined 
by the solution of the BBG equation [20], Vg eff must be 
calculated self-consistently with the G matrix and the 
s.p. potential Uk on the basis of BBG equation. It is clear 
from Eq. (|3"4")l that the effective force rising from the TBF 
in nuclear medium is density and temperature dependent 
through the defect function. A detailed description and 
justification of the method can be found in Ref. [20], in- 
cluding a discussion on the averaging procedure. The 
validity of such a procedure has been numerically tested 
in the comparison between the BHF EoS plus Z-diagrams 
with a meson exchange and the Dirac-BHF EoS, which 
are expected to be equal. The calculation [39| gives an 
impressive agreement between the two EoS's, although 
the TBF due to the Z-diagrams is averaged according to 

Eq. flMD- 

For V r 9 bare we adopt the Argonne Vis two-body inter- 
action |43|]. The TBF is constructed from the meson- 
exchange current approach [20] and contains virtual par- 
ticle (A and AT* (1440)) excitations and, in addition, rel- 
ativistic effects induced by the excitations of particle- 
antiparticle pairs. This description of the interaction is 
not completely consistent since, in principle, the two and 
three body forces should be derived from the same me- 
son parameters, but a recent calculation [30] replacing 
Argonne potential with Bonn potential [4J] and the TBF 
built up with the Bonn meson parameters substantially 
leads to the same results. 
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FIG. 7: Isotherms of symmetric nuclear matter (left side) and 
pure neutron matter (right side) as a function of the nucleon 
density at different temperature. 



B. Thermodynamics 

Let us start with symmetric nuclear matter. In the 
BHF approximation, the thermodynamic potential can 
be written 



9. 



fiSr + Wi, 



o'l = -T^ln[l + e-^-^/ T 



'N 



W L = -\Y.UT)U k {T), 



(35) 
(36) 

(37) 



where J7^j is the thermodynamic potential for a sys- 
tem of independent particle with single particle spec- 
trum e k = ^- + U k (T), and the W L the sum of all 
linked cluster diagrams to the lowest order in the hole 
line expansion [40j. Uk{T) is the selfconsistent mean field, 



u k( T ) = ^2{kk'\G\kk') A fk'(T). 



(38) 



The finite temperature BHF approximation suffers from 
the same difficulty as any strongly interacting Fermi sys- 
tem. The difficulty is the same as in the CDDM quark 
model, since in both cases the s.p. spectrum is density 
and temperature dependent. Whereas at T = the den- 
sity (or Fermi momentum) can be fixed, at T > the 
role of density is taken by chemical potential p (grand 
canonical ensemble), and one is forced to fix the chemi- 
cal potential at each iteration due to the presence of the 
Fermi distribution f k (T) = {1 + exp([e fe (T) - ^/T)}" 1 . 
Since this procedure does not converge [40], one should 
fix the density and invert the equation relating density 
and chemical potential, 



= T7E/*( T )- 



V 



dW L 



(39) 



which is not a viable task. The usual approximation is to 
drop out the derivative in the previous equation, which 
corresponds to the quasiparticle approximation above 
discussed within the CDDM quark model. So doing, the 
resulting /2 looses its meaning of chemical potential. In 
this approximations the energy and entropy densities are 
given by: 



E N 
Sn 



k 

-^E{^( T ) ln ^( T ) 



U k (T) 



(40) 



[l-fk(T)]\n[l-fk(T)]}. (41) 



After one calculates /2 in terms of p, the thermodynam- 
ics is developed from the free energy density Fj$ (p, T) = 
Ew{p, T) ~ TSn(p, T). The free energy per particle, cal- 
culated from BHF approximation, is depicted in Fig. [7] 



Due to the difficulty of extending the BHF code to very 
high temperature an extrapolation from the real numer- 
ical results to high T has been performed adopting the 
so called frozen approximation based on T-independent 
single particle spectrum, i.e. the latter is frozen at T = 0. 
This turns out to be a good approximation up to 10 ~ 20 
MeV H|. 

The relevant thermodynamical quantities, i.e. chemical 
potentials and pressure, are derived from free energy as 
follows 



1000 



P 



dR 



N 



dp 



P 



P 



2 d( F N /p) 
dp ■ 



(42) 



Let us consider asymmetric nuclear matter with 
baryon density p = p n + p p and asymmetry parameter 
P = (Pn — Pp)/p, where p n (p p ) is the neutron (proton) 
density. The baryon chemical potentials can be expressed 
as 



Pn = 



dp„ 



Pp 



T,p P 



dF N 
d Pp / a-.,„. 



(43) 



where En is the free energy density. 

Assuming the parabolic law for the latter, we get the 
simple expression for the chemical momentum isotopic 
shift 



Pn - Pp = 4 / 3E sym (p, T), 



(44) 



where E sym is the symmetry free energy density. The 
parabolic law is well satisfied at low density, but at high 
density additional terms of the /^-expansion must be con- 
sidered. 

In neutron star inner core nuclear matter is supposed 
to be in /^-equilibrium under the condition of charge neu- 
trality. Assuming that only electrons are present (the 
muon contribution is negligible), the two preceding con- 
ditions require 

Pe = Pn - Pp, (45) 

Qn = Pp - Pe- (46) 

Qn is the net charge density of nuclear matter. It is 
zero for pure neutral nuclear matter. For a given set of 
(p, Qn), w e can solve the chemical potentials /i n , p p , and 
fi e from the above equations. Then all other quantities 
can be obtained for a fixed temperature. In other words, 
all thermodynamic quantities can be regarded as a func- 
tion of the nucleon density p, charge density Qn, and 
temperature T. At zero temperature, for example, the 
energy density can be regarded as a function of p and 
Qn, i-e., En = E^(p, Qn)- With a similar approach as 
in the preceding section to obtain Eq. ([22]) . we can eas- 
ily show that the baryon chemical and charge chemical 
potentials of nuclear matter can be expressed as 



dE h 

dp„ 



Pn, 



5En 
5Q> 



= Pp- Pn- 



(47) 



The system turns out to be in a strongly isospin asym- 
metric state. The isotherms of free energy and pressure 
of nuclear matter in /3-equilibrium are shown in Fig. [5J 
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FIG. 8: The free energy per baryon (left y axis) and pressure 
(right y axis) of nuclear matter in /3-equilibrium. Three values 
of the temperature are considered. 



IV. PHASE DIAGRAM STRUCTURE AT ZERO 
AND FINITE TEMPERATURE 

Let us study the nuclear matter, consisting of nucleons 
and electrons, in equilibrium with a gas of u, d, s quarks 
and electrons. According to Glcndenning [4^, [47]], we 
assume the total charge conservation, in addition to total 
baryon and energy conservation. Now we first consider 
the case of zero temperature, and then extend to finite 
temperature. 

The conservation laws can be imposed by introducing 
the quark fraction \ defined as 



X = vjv. 



(48) 



where V is the total volume, V^ is the volume occupied 
by quarks. Then the total baryon density is 



pt = (i-x)p + x«, 

the total electric charge is 

Qt = (i-x)Qn +xQ q , 

and the total energy density is 

E t = (1 - X )E N + xE q , 



(49) 



(50) 



(51) 



where p, Qn, and En are, respectively, the baryon num- 
ber density, electric charge density, and energy density of 
nuclear matter, while n, Q q , and E q are the correspond- 
ing quantities of quark matter. Efa is a function of p and 
Qn, E q is a function of n and Q q , i.e., En = E^(p, Qn), 
E q = E q (n, Q q ). Differentiating Eq. (jSTj) . one obtains 



dE t . 



+X 



(1 - X)\ -r^—dp + T^-dQ N 



V dR 



dE 



dn 



dE t 



dQ. 



^dQ c 



<9Qn 

+ (Eq 



£ N )dx. (52) 
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On the other hand, differentiating Eqs. (I49[) and ([50)1 at 
a given pair of pt and Qt , we have 

(1 - x)dp = (p - n)&\ - Xdn, (53) 

(1 - x)dQ N = (On ~ Q q )dx - xdQ q . (54) 

To minimize Pt, we substitute Eqs. ([53]) and ((54|) into 
Eq. (52]). Then setting dP t = 0, we find 



dE N dE a dE h 



dp 



dE a 



dn ' <9Q N dQ q 



Pn = Pa 



where 



Pn = 



rr x. 9En ±n dEN 
-E N + p— \-Qn- 



E a 



dp 

dE q 
dn 



dQ> 
dE q 
l dQa 



(55) 

(56) 
(57) 



The conditions in Eq. ([55)) are nothing but the Gibbs 
ones, i.e., the baryon chemical potential, the charge 
chemical potential, and the pressure in nuclear and quark 
matter should be equal to each other to minimize the to- 
tal energy of the mixed phase. 

In the previous two sections, we have linked the 
baryon chemical potential and charge chemical potential 
to the respective constituent particle chemical potentials 
in Eqs. P^j) and (|4"7"|) . As application of these equalities, 
we immediately see that the first two equations in ([5f 
are equivalent to 



Pn = Pu + %Pd 



A^p — ^f^v 



Pd- 



p u = (2p p - Mn)/3, 

p d = (2p n - p p )/3. 



(58) 



In general, all other chemical potentials in quark sector 
can be related to p u and p d , e.g., p s = p d , p e = Pd~ Pu- 
Similarly, all chemical potentials in nuclear sector can 
be linked to p n and p p , e.g., p e — p n — p p . There- 
fore, Eq. (|58p means that we can choose either (p u ,pd) 
or (p n ,p p ) as the two independent chemical potentials. 
The latter can then be determined by solving the charge 
neutrality equation and the pressure balance equation for 
a given total baryon number or a given quark fraction. 

At finite temperature, we similarly have the phase 
equilibrium condition 



Pn 



p, 


(mechanical) , 


(59) 


Mq 


(chemical) , 


(60) 


^T 


(thcrmodynamical) . 


(61) 



Pn = 
Pn = T c 



The condition in Eq. (|61[) only tells us that the temper- 
ature in nuclear and quark sectors are equal, so we have 
a common temperature T. The chemical equilibrium 
condition in Eq. ([rJO]) is equivalent to that in Eq. (155)) . 
Therefore, we still have only two independent chemical 
potentials. For a given total density pt at a fixed tem- 
perature T, the two independent chemical potentials and 
the quark fraction \ can then be determined by solving 
the three equations in Eqs. ([4*5]) . (150)) with Q t = 0, and 

dFJD. 



Similar to the case at zero temperature, the lower crit- 
ical density p c \, which separates the nuclear and mixed 
phases, is defined by \ — 0, while the critical density p c i 
between the mixed and quark phases is determined by 
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FIG. 9: The free energy per baryon as a function of density in 
the nuclear (dashed curves), mixed (solid), and quark phase 
(dotted). The quark fraction is also shown on the right axis 
with dash-dotted lines. The curves with an open circle are 
at zero temperature while those with a full circle are for the 
temperature T = 20 MeV. 

In Fig. [HI we display, for the parameter \f~D = 170 
MeV, the density dependence of the energy per baryon 
in pure nuclear matter, in pure quark matter, and in the 
mixed phase. The quark fraction has also been depicted 
on the right y-axis. In Fig. [TUJ the corresponding pres- 
sure is reported. It is seen that the nuclear matter is 
the most favorite phase at lower densities, and the quark 
matter is the most stable phase at higher densities, while 
at intermediate densities, the mixed phase has the lowest 
energy. 

The quark baryon number density n and nuclear den- 
sity p are also plotted on the right axis of Fig. [TO] We 
see that the quark density is always higher than the nu- 
clear density. The transition from hadron phase to mixed 
phase occurs at the density a bit less than 0.15 fm , well 
below the saturation density. But it is hard to observe 
in terrestrial laboratories, since the nuclear matter so far 
realized in exotic nuclei or heavy-ion collisions is much 
less neutron rich. The transition from mixed phase to 
pure quark phase occurs at the total density p% — 0.85 
fm~ 3 , where the nuclear density is only p — 0.64 fm -3 . 
The density range of mixed phase is only slightly depend- 
ing on the temperature, at least in temperature range of 
interest in neutron stars. 

The critical densities depend on the parameter D. In 
the left panel of Fig. [Tl] both nuclear critical density 
(the solid line, separating the pure nuclear phase and the 
mixed phase) and the quark critical density (the dashed 
line, delimiting the pure quark phase) are displayed as a 
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FIG. 10: The pressure in nuclear (dashed curve), mixed 
(solid), and quark phase (dotted) vs density at temperature 
T — (cures with an open circle) and 20 MeV (cures with an 
solid circle), the nuclear density p and the quark density n in 
the mixed phase are also shown on the right axis. 



function of D 1 / 2 . If D 1 / 2 < 161.6 MeV, the two critical 
densities approach zero, and accordingly SQM is abso- 
lutely stable. When 161.6 MeV < D 1 / 2 < 162.5 MeV, 
mixed phase can exist at any lower densities. Only when 
D 1 / 2 > 162.5 MeV, nuclear matter is more stable at lower 
densities. If we have only two flavor quarks in the quark 
sector, the critical densities are usually higher. In the 
same panel we also plot the lower critical density for the 
two-flavor case. Because we know in our real world the 
two-flavor quark matter does not exist below the satura- 
tion density, D 1 / 2 should be on the right of the first full 
dot at £> 1/2 w 168 MeV (the intersection of the dotted 
and dot-dashed lines) in Fig. [TTJ On the other hand, to 
let SQM have a chance to appear below the saturation 
density, D 1 ' 2 should be on the left of the second full dot 
where D 1 / 2 = 171.3 MeV. In plotting Fig. M and E01 we 
adopt D 1 / 2 = 170 MeV. The temperature dependence of 
the lower and higher critical densities is also plotted in 
Fig. [TT] (right panel). The two lines at fixed D mark the 
boundaries of the three phases. 



V. PROPERTIES OF HYBRID STARS 

With the equation of state that has the mixed and/or 
quark phase derived in Sec. IIV1 we are ready to study 
the structure of hybrid stars by solving the Tolman- 
Oppenheimer-Volkov equation 



dP _ GmE (1 + P/E)(l + 47rr 3 P/m) 
dr r 2 1 — 2Gm/r 



(62) 



where G = 6.707 x 10~ 45 MeV" 2 is the gravitational 
constant, r is the distance from the center of the star, 
E = E(r) and P = P(r) are the energy density and 
pressure at the radius r, respectively. The subsidiary 
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FIG. 11: Left panel: The critical density of nuclear matter 
to quark matter as a function of the confinement parameter 
D. The horizontal line is the nuclear saturation density. The 
dotted line is the quark critical density in the case of two- 
flavor quark system. Right panel: Phase diagram for three 
values of the confinement parameter and m s o = 95 MeV. 



condition is 



dm/ 'dr = 4irr E 



(63) 



with m = 771 (r) being the mass within the radius r. 

At variance with pure nuclear or quark stars, a hybrid 
star contains pure quark matter in the core, pure nuclear 
matter near the outer part, and, in between, a mixed 
phase of the quark and nuclear matter. In this case, 
therefore, we must use the EoS in whole density range. 




nuclear matter 
D 1 ' 2 = 170MeV 
D 1 ' 2 = 180MeV 
D 1 ' 2 = 190MeV 
D 1 ' 2 = 200MeV 



10 12 14 1 
RADIUS (km) 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 
CENTRAL DENSITY (fm 3 ) 



FIG. 12: Mass-radius relation of hybrid stars at m s o — 95 
MeV for four values of the confinement parameter D 1 ' 2 of 
170, 180, 190, 200 MeV . 

The resulting gravitational mass for the hybrid star is 
plotted in Fig. [T2J as a function of both radius and cen- 
tral density, for four values of the confinement parameter 
D 1 ' 2 of 170, 180, 190, 200 MeV. The main effect of the 
phase transition is, as expected, a large reduction of the 
maximum mass due to softening of the EoS. Concern- 
ing the confinement parameter D, we observe a slight 
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decrease of the maximum mass when the D value be- 
comes smaller, from 1.74M Q at D 1 / 2 = 200 MeV down 
to 1.6, 1.55, 1.52 M Q corresponding to D 1 / 2 = 190, 180, 
170 MeV. This can be easily understood in this way: 
since quark phase is occurring at lower density for smaller 
values of D, sometimes even less than the nuclear satu- 
ration density (for example only 0.15 fm -3 for Z? 1 / 2 = 
170 MeV), then the quark population of the star would 
be more numerous, since the stronger softening of the 
EoS may only support less gravitational mass. In par- 
ticular, the most favorable case of D 1 ^ 2 = 180 and 190 
MeV, when the quark confinement appears at around 2po 
consistent with the heavy ion experiment, the predicted 
maximum mass can match very well the S branch of neu- 
tron stars mentioned by Haensel et al. with an error less 
than 5% 0. 
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FIG. 13: Temperature effect on the mass-radius relation of 
hybrid stars. 

In order to study the temperature effect on the mass 
of hybrid stars, we plot the star mass, vs the star ra- 
dius or central density, at T — and T — 20 MeV in 
Fig. [13J It is seen that the temperature influence on the 
maximum mass is very limited. Otherwise the effect is a 
quite strong increase of the NS radius for a fixed amount 
of gravitational mass. But, larger is the mass smaller is 
the radius variation, similarly to a few other calculations 

El. 



A. Kaon condensation 

The condensation of K~ mesons in neutron stars is 
widely discussed in the literature (see Refs. [49|, [5(| and 
references therein quoted). In dense matter, the conden- 
sation of K~ mesons is originated by the reaction 



K~ 



(64) 



if the effective mass of the K~ drops below the chemi- 
cal potential of the electron, this reaction becomes possi- 
ble in dense matter, indicating the presence of kaon con- 
densation. Since almost all the studies of kaon medium 



properties [5l| suggest the consistent picture that the at- 
traction from nuclear matter would bring the K~ mass 
down, so the threshold condition for the onset of K~ 
condensation [i e = m* Kl which follows from Eq. (I64p . 
could be fulfilled in the center of neutron stars [5l| at 
p It 3po- However, the deconfinement phase transition 
from the hadronic phase to the quark phase occurs also 
at rather low density, which leads to the onset of mixed 
phase in neutron stars. Then it is very interesting to ex- 
plore how the quark deconfinement affects the K _ con- 
densation threshold. 
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FIG. 14: The effective kaon mass and the electron chemi- 
cal potential are shown in normal nuclear matter and also in 
mixed phase at m 3 o = 95 MeV for three selected values of the 
confinement parameter D 1/2 of 180, 200, 225 MeV. 

We take the antikaon dispersion relation constrained 
by the heavy-ion data as empirical indication of an at- 
tractive antikaon potential in dense matter [5l| , and com- 
bine this with the BHF model of nuclear matter together 
with the above quark model without bag constant, to cal- 
culate the effective kaon mass in hybrid stars. The result 
is illustrated in Fig. [Mj where the effective kaon mass 
and the electron chemical potential are shown in normal 
nuclear matter and especially in mixed phase for three 
selected values of the confinement parameter D. One 
sees that in normal nuclear matter, K~ medium mass 
decreases with increasing density and meets the electron 
chemical potential at ~ 0.6 fm -3 (the solid bullet in Fig. 
14), so the kaon condensation would be present for sure 
in this case; However, once the quark phase sets in at 
the total density of 0.26, 0.4, 0.53 fm~ 3 respectively for 
D 1 / 2 = 180, 200, 225 MeV, the decreasing speed of the 
kaon effective mass in the matter slows down a little bit. 
In addition, more conclusively, the electron population 
begin to decrease instead, for the reason that the elec- 
tric charge neutrality can be achieved more efficiently 
through the charged quarks themselves. As a result the 
threshold condition of K~ condensation is much more 
difficult to satisfy, unless the confinement parameter D 
is chosen to be extremely high (at least D 1 / 2 = 225 MeV) 
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when the presence of quark phase is pushed to very high 
densities, then one may expect that m* K can finally equal 
to /i e (which is indicated with a circle). Those high val- 
ues of D may even not be realistic, therefore we would 
like to conclude that the inclusion of quark phase may 
make the kaon condensation impossible in neutron stars, 
or at least hinder it very strongly. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper the neutron star inner structure has been 
studied. The lack of strong observational constraints de- 
mands for sophisticated models of the NS composition 
and interaction mechanisms. In this study we only in- 
cluded baryons and quarks in equilibrium with leptons 
and kaons. Since in the NS interior high baryon densi- 
ties are reached, NN and nucleonic excitations are ex- 
pected to play a major role. Their effects can be in- 
corporated in a three body force. The baryon EoS in 
weak coupling equilibrium with electrons was derived 
within the BBG theory suitably extended so to include 
the three body force. In the quark sector, we adopted the 
semiphenomenological CDDM quark model, which ex- 
hibits a confinement mechanism alternative to the crude 
MIT bag model. Furthermore, in contrast with the ex- 
tension of the MIT model, where the density dependence 
is introduced artificially [52| , CDDM quark model shapes 
its densit y dep endence in agreement with chiral require- 
ments QJ, [j]f. 

Hadrons and quarks in /3-equilibrium were considered 
at zero and finite temperature. Since typical tempera- 
tures of protoneutron stars are as high as 40 MeV and 
beyond, rigorous thermodynamic potentials have to be 
derived. This was possible with CDDM quark model, 
due to its simplicity, whereas some approximations are 
needed to calculate the hadron phase in the Brueckncr 
theory. The latter is still a main drawback of the finite 
temperature microscopic theory of strongly interacting 
Fermi systems. 

The transition from the low density hadron phase to 
high density quark phase in beta equilibrium was studied 
under the Glendenning hypothesis of total charge neu- 
trality. The Gibbs construction enabled to follow the 
evolution of the mixed hadron-to-quark phase, varying 
the temperature T and the confinement parameter D. 
The EoS, in terms of pressure, energy density and chem- 
ical composition, was constructed as a function of T and 
D. Moreover the phase diagram T-D was also depicted. 
The transition density from hadron to mixed phase is 
strongly dependent on the confinement parameter and, 
for D small enough, it becomes lower than the nuclear 
saturation density. However, at that point nuclear mat- 
ter is so strongly neutron rich that it hard to imagine 
that it can be observed in terrestrial laboratory exper- 
iments. At high temperature the transition density is 
comparatively smaller. 

The TOV equations have been solved in the above dis- 



cussed model of neutron stars. The transition to quark 
matter produces a strong reduction of the maximum 
mass, from M = 2.3M Q to 1.5M Q for the lowest value 
of the confinement parameter D, but the corresponding 
radius is not changed, being in both cases about 10 km. 
The other M-R configurations ly in a range R ~ 10 ~ 13 
km depending on the value of D. With increasing D, the 
maximum mass increases from 1.5 to 1.7 solar masses. 
Other quark models have been adopted to describe the 
NS transition to a deconfined phase (see, for instance, 
Refs. [Hill SI). Except for the Nambu-Jona Lasinio 
model which exhibits instability of neutron stars [24| , the 
other models point to a softening of the EoS of nuclear 
matter. The MIT model with density-dependent bag 
constant predicts a beta-stable quark phase quite similar 
to our model [4g], but the comparison on the NS struc- 
ture turns out to be difficult, since in that calculation the 
transition to quark phase is built on top of hyperonized 
nuclear matter. 

The effect of temperature on the maximum mass is 
negligible up to T — 20 MeV, but as expected the heated 
system can support much less gravitational mass in the 
same volume for all M-R configurations. Thus, an evo- 
lutionary study of isolated neutron stars would show a 
strong compression from the newborn phase to the long 
era phase as fast as the star cools down from 40 MeV to 
approximately zero. 

Beyond the confinement parameter, the transition to 
quark phase is interrelated to other properties, in par- 
ticular the possible onset of kaon condensation. In our 
analysis we concluded that the quark phase could not be 
compatible with the kaon condensation in neutron stars, 
or at least could hinder it very strongly. 

The model developed in this paper misses some im- 
portant aspects. The first one is the effect of neutrinos. 
Neutrinos play the major role in new born neutron stars, 
when they are still trapped in the interior. In this case 
the neutronization process is still hindered and the EoS 
is that of symmetric nuclear matter. Moreover, the on- 
set of kaon condensation is shifted to higher densities [53| 
and probably kaons have no more any chance to compete 
with quarks. 

The second one is the inclusion of hyperons and their 
competition with other mechanisms such as kaon con- 
densation. In Refs. [48] the hybrid stars are studied by 
including hyperons in the hadron phase, which makes the 
hadron EoS very soft and brings the maximum mass to 
M = 1.25Mq. The introduction of the quark phase rises 
up again the maximum mass to 1.5Mq, a value which 
is consistent with our prediction without hyperons. On 
the other hand, the appearance of hyperons much de- 
pends on the threshold of kaon production. Therefore 
the interplay between hyperons and kaons turns out to 
be quite important and deserves additional investigation, 
as soon as more reliable empirical inputs will be available, 
especially on the hyperon-nucleon and hyperon-hyperon 
interaction. 
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APPENDIX A: THERMODYNAMICS WITH 
CONFINEMENT BY THE DENSITY 
DEPENDENCE OF QUARK MASSES 

Let's start from the fundamental thermodynamic 
differentiation relation d{VE) = Td(VS) - PdV + 
^nPidiVrii), where S is the entropy density, and 
(anti)particles are assumed to be uniformly distributed 
in a volume V. By using the free energy density 
F = E - TS, it becomes d(VF) = -VSdT - PdV + 
J2, md(Vm), or, equivalent^, dF = -SdT + (-P -F + 
J2i ^i n i)<iV/V +J2i Mjdn.i. Because of the uniformity, the 
free energy density has nothing to do with the volume. 
We thus have 



P = -F + ^fljTli, 
i 

dF = -SdT + ^Hidm 



(Al) 
(A2) 



At finite temperature, we should consider both par- 
ticles and anti-particles, and particles/anti-particles are 
not always located bellow the Fermi energy. Therefore, 
The net particle number densities can be expressed as 



Hi = Ui(T, p*,mi) = nf - n\ 



(A3) 



where the superscripts + indicates particles, and signi- 
fies antiparticles: 



nf = gi 



1 



p 2 dp 



1 _|_ e (vV+™?TM*)/T 2tt 2 



(A4) 



In Eqs. (|A4[) and (|A3|) . p* are effective chemical po- 
tentials of respective particles. If quark densities are not 
density and/or temperature dependent, p* are nothing 
but the actual chemical potentials. In our present case, 
however, quark masses depend on both density and tem- 
perature, in order to include the strong interaction be- 
tween quarks. Therefore, the real chemical potentials 
should be derived according to fundamental thermody- 
namic laws. 

Eq. (|A3[) gives, implicitly, p* as a function of T, rn, 
and rrii, i.e., 



Mi 



H*(T,ni,mi). 



(A5) 



To determine the thermodynamic properties, we need 
to give a characteristic function. At zero temperature, 
we use the energy density in Eq. @ due to zero entropy. 
Now the temperature T and the densities rii are chosen as 
the independent system variables, we should, therefore, 
choose the free energy in Eq. (|29|) as the characteristic 
function. 

Differentiation of Eq. (f2"9")) gives 



dF = £ 



m ,dFd K \ dT+ mo^ 



dT dp* dT 
8F dF dp* 



dp* dm 



drrii dp* dm 

drnj irp I y-v drrij 



drrij 



(A6) 



Applying drrii — ^y-dT+Y\. ¥^drij, then comparing 
the corresponding expression with Eq. (IA2|) we immedi- 
ately have 



(M 



dF dp* 
dp,* dn t 

and 

* = -E 



E 



dFj dF 3 dp* \ dm, 



drrij dp* drrij I drii 



dF dFdpl (m dFdti 
dT dp* dT \drrn dp* dm., 



(A7) 



drrij 



dT 
(A8) 



To simplify the expressions, we differentiate Eqs. (IA3 
and llAlfll to get 



A drli Ar 

1 = ~dT 



drij 

dpi 



dp* 



* dp* dp* 

which implies 
dp* drii drii dp* drii 



drij 
dm, 



drrii, 
dm,:, 



drrii 

dp* drij 

dT dp* dT ' dn t dp* ' drm dp* 

Defining 
Q = J2^o,i(T,p*,mi) = J2[ n ' 



(A9) 
(A10) 

dm 



+ 

0,i 



n 



0,/J 



drrii 
(All) 

(A12) 



with 



a 



0,! 



g{T 

2tt 2 



In 



1 + p-iV^+^TuD/T 



P dp, 



. (A13) 
then we can write F^ = fl oi + p*rii. Substituting this 

into Eq. (TA8J) and (|A7|) . then applying Eq. (jAllj) and 

rii — ~dflo t i/dp*, we have 



Mi = M* + E 



<9£7o drrij 
drrij drii 



and 



dflo sr-^ dflo drrii 



dT 



dm, dT 



(A14) 



(A15) 



Obviously, the free energy density can be given as 

dn 



F = n -J2 



Mt 



dn\ 



(A16) 



The energy density is obtained by E = F + TS, giving 
E = VL -^ l x i —-T—-T^ — — . (A17) 



dT 



dm, dT 



And the pressure is obtained by substituting Eq. (|A14[) 
into Eq. ([AT]) : 



p=-n +X/' 



dflg drrij 
' drrij drii 



(A18) 



Eqs. (|A15|) - (|A18|) here are in complete accordance with 
the Eqs. (58)-(61) in [19[ if one regards the rii there as 
\x\ and £1 as Slo which were not explicitly stated. The 
expression in Eq. (|A14|) is special in the present paper 
because we need the real chemical potential to investigate 
the mixed phase. 

The real thermodynamic potential density of the sys- 
tem is 



tt = F -^ / A ' n » = ^0 ~ z2 '' 



1,3 



dH, drrij 
drrij drii 



(A19) 
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In the above derivation, we choose volume V, the 
temperatureT, and the particle number densities n, as 
the independent system variables. In this case, the free 
energy is the characteristic function from which we get 
the complete set of thermodynamic functions. For this 
purpose we have defined the intermediate variables /i* 
in Eqs. (|A4|) and ([50)1 . Because the quark matter we 
are considering is a strongly interacting system, the rela- 
tions between the chemical potentials and the densities 
are, in principle, not the same as those of a free Fermi 
gas. However, with the 'effective' chemical potentials fi*, 
the densities and the free energy are really of the same 
form as those of a non-interacting Fermi gas. This is 
what the 'effective' means. The actual chemical poten- 
tials of each type of particles are determined from the 
fundamental thermodynamic equality (|A2|) which results 
in Eq. (|A14|) . 
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